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Abstract 

In this paper, we focus on the problem of existence and computing of small and large 
stable models. We show that for every fixed integer k, there is a linear-time algorithm to 
decide the problem LSM (large stable models problem): does a logic program P have a 
stable model of size at least \P\ — fe? In contrast, we show that the problem SSM (small 
stable models problem) to decide whether a logic program P has a stable model of size 
at most k is much harder. We present two algorithms for this problem but their running 
time is given by polynomials of order depending on k. We show that the problem SSM is 
fixed-parameter intractable by demonstrating that it is VT[2]-hard. This result implies that 
it is unlikely an algorithm exists to compute stable models of size at most k that would 
run in time 0(m c ), where m is the size of the program and c is a constant independent 
of k. We also provide an upper bound on the fixed-parameter complexity of the problem 
SSM by showing that it belongs to the class W[3]. 



1 Introduction 



The stable model semantics by Gelfond and Lifschitz ( Gelfond fe Lifschitz, 198cf ) is 
one of the two most widely studied semantics for normal logic programs, the other 



one being the well-founded semantics by Van Gelder, Ross and Schlipf (Van Gelder 



et ai, 1991). Among 2-valued semantics, the stable model semantics is commonly 



regarded as the one providing the correct meaning to the negation operator in logic 
programming. It coincides with the least model semantics on the class of Horn pro- 
grams, and with the well-founded semantics and the perfect model semantics on 
the class of stratified programs (Apt et aL, 1988). In addition, the stable model 



semantics is closely related to the notion of a default extension by Reiter (Marek 
|fe Truszczyhski, 1989| [Bidoit fe Froidevaux, 199l"| ). Lo gic programming with stable 
model semantics has applications in knowledge representation, planning and reason- 
ing about action. It was also recently proposed as a computational paradigm well 
suited for solving combinatorial optimization and constraint satisfaction problems 



( [Marek fe Truszczyhski, 1999| ; [Niemela, 1999Q . 

Before we proceed, we will recall the definition of a stable model of a logic pro- 



gram, and some related terminology and properties. The reader is referred to (Marek 
|fe Truszczyhski, 1993 ) for a more detailed treatment of the subject. In the paper 



* This is a full version of an extended abstract presented at the International Conference on Logic 
Programming, ICLP-99 and included in the proceedings published by MIT Press. 
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we deal only with the prepositional case. A logic program rule is an expression r of 
the form 

r = a <— ,6 s ,not(ci), . . . ,not(c t ), 

where a, biS and CiS are propositional atoms. The atom a is called the head of r 
and is denoted by h(r). Atoms 6, and c, form the body of r. The set {bi,...,b s } is 
called the positive body of r (denoted by b + (r)) and the set {ci, . . . , c t } is called 
the negative body of r (denoted by 6 — (r)). A Zogic program is a collection of rules. 
For a logic program P, by At{P) we denote the set of atoms occurring in its rules 
and by h(P) — the set of atoms appearing as the heads of rules in P. We will also 
denote the size of P, that is, the total number of occurrences of atoms in P, by 
size(P). Throughout the paper we use n to denote the number of atoms in a logic 
program P, and m to denote the size of P. 

A set of atoms M C At(P) satisfies a rule r if h(r) e M, or if b+(r) \ M 7^ 0, 
or if 6 _ (r) PI M ^ 0. A set of atoms M C ^4i(P) is a model of a program P if M 
satisfies all rules of P. 

A logic program rule r is called Horn if b~ (r) = 0. A Horn program is a program 
whose every rule is a Horn rule. The intersection of two models of a Horn program 
P is a model of P. Since the set of all atoms is a model of P, it follows that 
every Horn program P has a unique least model. We will denote this model by 
LM(P). The least model of a Horn program P can be constructed by means of the 
van Emden-Kowalski operator Tp ( van Emdcn fc Kowalski, 1976Q . Given a Horn 



program P and a set of atoms M C P, we define 

Tp(M) = {a: a 6 X , . . . , b s e P, and {61, ...,&,} C M}. 
We also define 

T P (M) = 0, and T p +1 (M) = T P (T P (M)). 

Since the operator Tp is monotone, the sequence Tp(0) is monotone and its union 
yields the least model of a Horn program P. That is, 

00 

LM(P) = {jTp($). 
t=o 

If P is finite, the sequence stabilizes after finitely many steps. 

For a logic program rule r, by horn(r) we denote the rule obtained from r by 
eliminating all negated atoms from the body of r. If P is a logic program, we define 
horn(P) = {horn(r):r 6 P}. 

Let P be a logic program (possibly with rules containing negated atoms). For a 
set of atoms M C At(P) we define the reduct of P with respect to M to be the 
program obtained by eliminating from P each rule r such that b~ (r) n M 7^ (we 
call such rules blocked by M), and by removing negated atoms from all other rules 
in P. The resulting program is a Horn program. We will denote it by P M . As a 
Horn program, P M has the least model LM(P M ). If M = LM(P M ), M is a stable 
model of P. Clearly, if M is a stable model of P, M C h(P). Both the notion of the 



reduct and of a stable model are due to Gelfond and Lifschitz (Gelfond & Lifschitz 
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In the paper we restrict our attention to programs whose rules do not contain 
multiple positive occurrences of the same atom nor multiple negative occurrences of 
the same atom in the body. It is clear that adopting this assumption does not limit 
the generality of our considerations. Repetitive occurrences can be eliminated in 
linear time (in the size of the program) and doing so does not affect stable models 
of the program. 

If M is a stable model of P, each rule r such that b + (r) C M and b~(r) D M = 
(that is, such that M satisfies its body), is called a generating rule for M. Clearly, 
if M is a stable model of P, it is also a stable model of the program consisting of 
all rules in P that are generating for M. 

There are several ways to look at the search space of possible stable models of 
a program P. The most direct way is to look for stable models by considering all 
candidate subsets of h{P). For each candidate subset M C h(P), one can compute 
the corresponding reduct P M , its least model LM(P M ), and check the equality 
M = LM(P ) to decide whether M is stable. An alternative way is to observe 
that stable models are determined by subsets of the set of atoms appearing negated 
in P. Indeed, let us denote this set by Neg(P) and let us consider sets M C At(P) 
and B C Neg(P). Let B' = Neg(P)\B. Then, M is a stable model of P if and only 
if M = LM(P B '), BC\M = and B' C M. Thus, the existence of stable models can 
be decided by considering subsets of Neg(P). Finally, one can consider the search 
space of all subsets of P itself, and regard each such subset as a candidate for the 
set of generating rules of a stable model. Indeed, if M C At(P) and P' C P, then 
M is a stable model of P if and only if M = h(P'), P' is the set of all generating 
rules for M in P and M = LM(horn(P')). 

The problem with the stable model semantics is that, even in the propositional 
case, reasoning with logic programs under the stable model semantics is compu- 
tationally hard. It is well-known that deciding whether a finite propositional logic 
program has a stable model is NP-complete ( Marek fc Truszczyhski, 1991 ). Conse- 
quently, it is not at all clear that logic programming with the stable model semantics 
can serve as a practical computational tool. 

This issue can be resolved by implementing systems computing stable models and 
by experimentally studying the performance of these systems. Several such projects 
are now under way. Niemela and Simons (Niemela & Simons, 1996) developed a 
system, smodels, for computing stable models of finite function symbol-free logic 
programs and reported very promising performance results. For some classes of 
programs, smodels decides the existence of a stable model in a matter of seconds 
even if an input program consists of tens of thousands of clauses. Encouraging 
results on using smodels to solve planning problems are reported in (Niemela, 1999). 
Another well- advanced system is DeReS (Cholewihski et al., 1996), designed to 
compute extensions of arbitrary propositional default theories but being especially 
effective for default theories encoding propositional logic programs. Finally, systems 



capable of reasoning with disjunctive logic programs were described in (Eitcr et al. 



1997) and (Aravindan et al., 1997). 



However, faster implementations will ultimately depend on better understanding 
of the algorithmic aspects of reasoning with logic programs under the stable model 
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semantics. In this paper, we investigate the complexity of deciding whether a finite 
prepositional logic program has stable models of some restricted sizes. Specifically, 
we study the following two problems (|P| stands for the number of rules in a logic 
program P): 

LSM (Large stable models) Given a finite propositional logic program P and an 
integer k, decide whether there is a stable model of P of size at least |P| — k. 

SSM (Small stable models) Given a finite propositional logic program P and an 
integer k, decide whether there is a stable model of P of size no more than k. 

Inputs to the problems LSM and SSM are pairs (P, k) , where P is a finite propo- 
sitional logic program and k is a non-negative integer. Problems of this type are 
referred to as parametrized decision problems. By fixing a parameter, a parame- 
terized decision problem gives rise to its fixed-parameter version. In the case of 
problems LSM and SSM, by fixing k we obtain the following two fixed-parameter 
problems (k is now no longer a part of input): 

LSM(k) Given a finite propositional logic program P, decide whether P has a stable 

model of size at least |P| — k. 
SSM(k) Given a finite propositional logic program P, decide whether P has a stable 

model of size at most k. 

The problems LSM and SSM are NP-complete. It follows directly from the NP- 



completeness of the problem of existence of stable models (Marek & Truszczynski, 



1991). But fixing k makes a difference! Clearly, the fixed-parameter problems SSM(k) 
and LSM(k) can be solved in polynomial time (unlike the problems SSM and LSM 
which, most likely, cannot). Indeed, consider a finite propositional logic program P. 
Then, there are 0(n k ) subsets of At(P) (in fact, as pointed out earlier, it is enough 
to consider subsets of h(P) or Neg(P)) of cardinality at most k (we recall that in 
the paper n stands for the number of atoms in P). For each such subset M, it can 
be checked in time linear in m — the size of P — whether M is a stable model of 
P. Thus, one can decide whether P has a stable model of size at most k in time 
0(mn k ). 

Similarly, there are only 0(|P| fc ) subsets of P of size at least |P| — k. Each such 
subset is a candidate for the set of generating rules of a stable model of size at least 
P| — k (and smaller subsets, clearly, are not). Given such a subset R, one can check 
in time 0(m) whether R generates a stable model for P. Thus, it follows that there 
is an algorithm that decides in time 0(m\P\ k ) whether a logic program P has a 
stable model of size at least |P| — k. 

While both algorithms are polynomial in the size of the program, their asymptotic 
complexity is expressed by the product of the size of a program and a polynomial 
of order k in the number of atoms of the program or in the number of rules of 
the program. Even for small values of k, say for k > 4, the functions mn k and 
TO |p|fc g row very fast with m = size(P), n — \At(P)\ and |P|, and render the 
corresponding algorithms infeasible. 

An important question is whether algorithms for problems SSM{k) and LSM(k) 
exist whose order is significantly lower than k, preferably, a constant independent of 
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fc. The study of this question is the main goal of our paper. A general framework for 
such investigations was proposed by Downey and Fellows (Downey & Fellows, 1997). 
They introduced the concepts of fixed-parameter tractability and fixed-parameter 
intractability that are defined in terms of a certain hierarchy of complexity classes 
known as the W hierarchy. 

In the paper, we show that the problem LSM is fixed-parameter tractable and 
demonstrate an algorithm that for every fixed fc decides the problem LSM(k) in lin- 
ear time — a significant improvement over the straightforward algorithm presented 
earlier. 

On the other hand, we demonstrate that the problem SSM is much harder. We 
present an algorithm to decide the problems SSM(k), for fc > 1, that is asymp- 
totically faster than the simple algorithm described above but the improvement is 
rather insignificant. Our algorithm runs in time 0(mn k ~ 1 ), an improvement only by 
the factor of n. The difficulty in finding a substantially better algorithm is not coin- 
cidental. We provide evidence that the problem SSM is fixed-parameter intractable. 
This result implies it is unlikely that there is an algorithm to decide the problems 
SSM(k) whose running time would be given by a polynomial of order independent 
of k. 

The study of fixed-parameter tractability of problems occurring in the area of 
nonmonotonic reasoning is a relatively new research topic. Another paper that pur- 
sues this direction is ( Gottlob et al, 1999| ). The authors focus there on parameters 
describing structural properties of programs and show that in some cases, fixing 
these parameters leads to polynomial algorithms. 

Our paper is organized as follows. In Section @, we recall basic concepts of the 



theory of fixed-parameter intractability by Downey and Fellows (Downey & Fellows 



1997 ). The following two sections present the algorithms to decide the problems 



LSM and SSM, respectively. The next section focuses on the issue of fixed-parameter 
intractability of the problem SSM and contains the two main results of the paper. 
The last section contains conclusions and open problems. 



2 Fixed-parameter intractability 

This section recalls basic ideas of the work of Downey and Fellows on fixed-parameter 



intractability. The reader is referred to ( Downey fc Fellows, 1997 ) for a detailed 
treatment of this subject. 

Informally, a parametrized decision problem is a decision problem whose inputs 
are pairs of items, one of which is referred to as a parameter. The graph colorability 
problem is an example of a parametrized problem. The inputs are pairs (G, fc), where 
G is an undirected graph and fc is a non-negative integer. The problem is to decide 
whether G can be colored with at most fc colors. Another example is the vertex 
cover problem in a graph. Again, the inputs are graph-integer pairs (G, fc) and the 
question is whether G has a vertex cover of cardinality fc or less. The problems 
SSM and LSM are also examples of parametrized decision problems. Formally, a 
parametrized decision problem is a set L C S* x E*, where S is a fixed alphabet. 

By selecting a concrete value a £ E* of the parameter, a parametrized decision 
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problem L gives rise to an associated fixed-parameter problem L a — {x : (x, a) G 
L}. For instance, by fixing the value of k to 3, we get a fixed-parameter version 
of the colorability problem, known as 3-colorability. Inputs to the 3-colorability 
problem are graphs and the question is to decide whether an input graph can be 
colored with 3 colors. Clearly, the problems SSM(k) (LSM(k), respectively) are 
fixed-parameter versions of the problem SSM (LSM, respectively). 

The interest in the fixed-parameter problems stems from the fact that they are 
often computationally easier than the corresponding parametrized problems. For 
instance, the problems SSM and LSM are NP-complete yet, as we saw earlier, 
their parametrized versions SSM(k) and LSM(k) can be solved in polynomial time. 
Similarly, the vertex cover problem is NP-complete but its fixed-parameter versions 
are in the class P. To see this, observe that to decide whether a graph has a vertex 
cover of size at most fc, where A: is a fixed value and not a part of an input, it is 
enough to generate all subsets with at most k elements of the vertex set of a graph, 
and then check if any of them is a vertex cover. A word of caution is in order here. 
It is not always the case that fixed-parameter problems are easier. For instance, the 
3-colorability problem is still NP-complete. 

As we already pointed out, the fact that a problem admits a polynomial-time 
solution does not necessarily mean that practical algorithms to solve it exist. An 
algorithm that runs in time 0(iV 15 ), where N is the size of the input, is hardly 
more practical than an algorithm with an exponential running time (and may even 
be a worse choice in practice). The algorithms we presented so far to argue that 
the problems SSM(k), LSM(k) and the fixed-parameter versions of the vertex cover 
problem are in P rely on searching through the space of N k possible solutions 
(where N is the number of atoms of a program, the number of rules of a program, 
or the number of vertices in a graph, respectively). Thus, these algorithms are not 
practical, except for the very smallest values of k. The key question is how fast 
those polynomial-time solvable fixed-parameter problems can really be solved. Or, 
in other words, can one significantly improve over the brute-force approach? 

A technique to deal with such questions is provided by the fixed-parameter in- 



tractability theory of Downey and Fellows ( Downey fc Fellows, 1997 ). A parametrized 



problem L C E* x E* is fixed-parameter tractable if there exist a constant p, an 
integer function / and an algorithm A such that A determines whether (x, y) G L 
in time /(|j/|)|£| p (\z\ stands for the length of a string z G £*). The class of fixed- 
parameter tractable problems will be denoted by FPT. Clearly, if a parametrized 
problem L is in FPT, each of the associated fixed-parameter problems L y is solvable 
in polynomial time by an algorithm whose exponent does not depend on the value 



of the parameter y. It is known (see (Downey & Fellows, 1997)) that the vertex 
cover problem is in FPT. 

There is substantial evidence to support a conjecture that some parametrized 
problems whose fixed-parameter versions are in P are not fixed-parameter tractable. 
To study and compare complexity of parametrized problems Downey and Fellows 
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proposed the following notion of reducibilityg. A parametrized problem L can be 
reduced to a parametrized problem L' if there exist a constant p, an integer function 
q and an algorithm A that to each instance (x, y) of L assigns an instance (x', y') 
of V such that 

1. x' depends upon x and y and y' depends upon y only, 

2. A runs in time 0(g(|?/|)|x| p ), 

3. (x, y) 6 L if and only if (x', y') 6 L'. 

Downey and Fellows also defined a hierarchy of complexity classes called the W 
hierarchy: 

FPT C W[l] C W[2] C W[3] ... (1) 

The classes W[t] can be described in terms of problems that are complete for 
them (a problem D is complete for a complexity class £ if D G £ and every problem 
in this class can be reduced to D). Let us call a boolean formula t-normalized if it 
is of the form of product-of-sums-of-products ... of literals, with t being the number 
of products-of, sums-of expressions in this definition. For example, 2-normalized 
formulas are products of sums of literals. Thus, the class of 2-normalized formu- 
las is precisely the class of CNF formulas. We define the weighted t-normalized 
satisfiability problem as: 

WS(t) Given a i-normalized formula (p, decide whether there is a model of ip with 
exactly k atoms (or, alternatively, decide whether there is a satisfying valuation 
for ip which assigns the logical value true to exactly k atoms) 

Downey and Fellows show that for t>2, the problems WS(t) are complete for the 
class W[t]. They also show that a restricted version of the problem WS(2); 

WS%{2) Given a 3CNF formula p and an integer k (parameter), decide whether 
there is a model of ip with exactly k atoms 

is complete for the class W[l]. Downey and Fellows conjecture that all the implica- 
tions in (Q) are proper^. In particular, they conjecture that problems in the classes 
W[t], with t > 1, are not fixed-parameter tractable. 

In the paper, we relate the problem SSM to the problems WS(2) and WS(3) to 
place the problem SSM in the W hierarchy, to obtain estimates of its complexity 
and to argue for its fixed-parameter intractability. 



3 Large stable models 

In this section we will show an algorithm for the parametrized problem LSM that 
runs in time 0(2 k+k m), where (P, k) is an input instance and, as in all other 



The definition given here is sufficient for the needs of this paper. To obtain structural theorems 

a Subtler definition is needed This topi 



: goes beyond the scope of the present paper. The reader 
is referred to ( Downey & Fellows, 1997 ) for more details. 

If true, this conjecture would imply that in the context of fixed-parameter tractability there is 
a difference between the complexity of weighted satisfiability for 3CNF and CNF formulas. 



8 



M. Truszczyriski 



places in the paper, m = size(P). This result implies that the problem LSM is 
fixed-parameter tractable and that there is an algorithm that for every fixed k 
solves the problem LSM(k) in linear-time. 

Given a logic program P, denote by P* the logic program obtained from P by 
eliminating from the bodies of the rules in P all literals not (a), where a is not the 
head of any rule from P. The following well-known result states the key property 
of the program P* . 

Lemma 3.1 

A set of atoms M is a stable model of a logic program P if and only if M is a stable 
model of P* . 



Lemma 3.1 implies that the problem LSM has a positive answer for (P, k) if and 
only if it has a positive answer for (P*,k). Moreover, it is easy to see that P* 
can be constructed from P in time linear in the size of P. Thus, when looking for 
algorithms to decide the problem LSM we may restrict our attention to programs 
P in which every atom appearing negated in the body of a rule appears also as the 
head of a rule (that is, to such programs P for which we have Neg(P) C h(P)). 

By P k let us denote the program consisting of those rules r in P for which 
|fr~( r )| < k. We have the following lemma. 

Lemma 3.2 

Let P be a logic program such that Neg(P) C h(P). Let M C At(P) be a set of 
atoms such that \M\ > |P| - k. Then: 

1. M is a stable model of P if and only if M is a stable model of P k 

2. if M is a stable model of P k , then P k has no more than k + k 2 different 
negated literals appearing in the bodies of its rules. 

Proof: (1) Consider a rule r e P \ P k . Then |&~(r)| > k + 1 and, consequently, 
&-(r) nM^l Indeed, if &~(r) n M = 0, then |M U 6"(r) = |M| + \b~{r)\ > |P|. 
Since Neg(P) C ft(P), &~(r) C fc,(P). In addition, (both if we assume that M 
is a stable model of P and if we assume that M is a stable model of P fc ), we 
have M C h(P). Thus, 6~(r) U M C h(P). Now observe that |P| > |ft(P)|. Thus, 
|M U&~~(r)| < |ft(P)| < |P|, a contradiction. 

Since for every rule r G P \ P k we have 6 - (r) nM/|, it follows that (P fe ) M = 
P M . Hence, M = LM(P M ) if and only if M = LM{{P k ) M ). Consequently, M is a 
stable model of P if and only if M is a stable model of P k . 

(2) Let P' be the set of rules from P k such that r G P' if and only if &~(r) flM = 
(the rules in P' contribute to the reduct (P k ) M ) and let P" be the set of the 
remaining rules in P k (these are the rules that are eliminated when the reduct 
(P fc ) M is computed). Since Neg(P) C h(P), for every rule r £ P, & _ (r) C /i(P). 
Thus, lj{ 6_ ( r ) :r G p '} ^ /i(P) \ M. Since M C ft(P) (as M is a stable model of 
P k ) and |P| > \h(P)\, we have | U{ fo_ ( ? ') :? ' G p '}l < fc - Further, since |P'| > |M| > 
|P| -fc > |P fe |- fc : ii; follows that |P"| < fc. Consequently, | \J{b~(r):r G P"}| < fc 2 . 
Hence, the second part of the assertion follows. □ 
Let us now consider the following algorithm for the problem LSM(k) (the input 
to this algorithm is a logic program P). 
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1. Eliminate from the input logic program P all literals not (a), where a is not 
the head of any rule from P. Denote the resulting program by Q. 

2. Compute the set of rules Q k consisting of those rules r in Q for which |6~(r)| < 
k. 

3. Decide whether Q k has a stable model M such that \M\ > \Q\ — k. 

This algorithm reports YES if and only if the program Q k has a stable model M 
such that \M\ > \Q\ — k. By Lemma |3~2|, that happens precisely if and only if Q has 
a stable model M such that \M\ > \Q\ — k. This last statement, by Lemma 3.1, is 



equivalent to the statement that P has a stable model M such that \M\ > \P\ — k. 
In other words, our algorithm correctly decides the problem LSM(k) . 

Let us notice that steps 1 and 2 can be implemented in time O(m), where the 
constant hidden by the "big O" notation does not depend on k. To implement 
step 3, let us recall that every stable model of a logic program is determined by 
some subset of the set of atoms that appear negated in the program (each such 
subset uniquely determines the reduct, as we stated in the introduction; see also 
( Bondarenko et ai, 1993 )). By Lemma 3.2, the set of such atoms in the program 
Q k has cardinality at most k + k 2 . Checking for each subset of this set whether it 
determines a stable model of Q k can be implemented in time 0(size(Q k )) — 0(m). 
Consequently, our algorithm runs in time 0(2 + m) (with the constant hidden by 
the "big O" notation independent of k). 

Theorem 3.3 

The problem LSM is fixed-parameter tractable. Moreover, for each fixed k there is 
a linear-time algorithm to decide whether a logic program P has a stable model of 
size at least \P\ — k. 



4 Computing stable models of size at most k 

In the introduction we pointed out that there is a straightforward algorithm to 
decide the problem SSM(k) that runs in time 0(mn k ), where m — size(P) and 
n = \At(P)\. For k > 1 (the assumption we adopt in this section), this algorithm 
can be slightly improved. Namely, we will now describe an algorithm for the problem 
SSM(k) that runs in time 0(F(k)mn k ~ 1 ), where F is some integer function. Thus, 
if k is fixed and not a part of the input, this improved algorithm runs in time 
0(mn fe_1 ). 

We present our algorithm under the assumption that input logic programs are 
proper. We say that a logic program rule r is proper if: 

(PI) h(r) i b+(r), and 
(P2) b+(r)nb-(r) = 

We say that a logic program P is proper if all its rules are proper. Rules that violate 
at least one of the conditions (PI) and (P2) (that is, rules that are not proper) have 
no influence on the collection of stable models of a program as we have the following 



well-known result (see, for instance, (Brass & Dix, 1997)). 
Lemma J^.l 
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A set of atoms M is a stable model of a logic program P if and only if M is a stable 
model of the subprogram of P consisting of all proper rules in P. 

It is easy to see that rules that violate (PI) or (P2) can be eliminated from a 
logic program P in time 0{m). Thus, the restriction to proper programs does not 
affect the generality of our discussion. 

For a proper logic program P and for a set A C At(P) of atoms, we define P(A) 
to be the program consisting of all those rules r of P that are not blocked by A (in 
other words, those that satisfy b~(r)nA — 0) and whose positive body is contained 
in A (in other words, such that b + (r) C A). 

Let P be a logic program and let A C At(P) be a set of atoms. A stable model 
M of P is called A-based if 

1. M is of the form A U {a}, where a 6 Ai(P) \ A, and 

2. M C LM(P{A) M ) (in other words, when computing LM(P M ), the derivation 
of A does not require that a be derived first). 

We have the following simple lemma. 

Lemma 4- -2 

Let k be an integer such that k > 1. A proper logic program P has a stable model 
of cardinality k if and only if for some A C At(P), with \A\ = k — 1, P has an 
A-based stable model. 



It follows from Lemma [4.2| that when deciding the existence of fc-element stable 
models, k > 1, it is enough to focus on the existence of A-based stable models. This 
is the approach we take here. In most general terms, our algorithm for the problem 
SSM(k) consists of generating all subsets A C At(P), with \A\ <k~ 1, and for each 
such subset A, of checking whether P has an A-based stable model. This latter task 
is the key. 

We will now describe an algorithm that, given a logic program P and a set 
A C At(P), decides whether P has an A-based stable model. To this end, we define 
P'(A) to be the program consisting of all those rules r of P such that: 

1. &"(r)ni = (r is not blocked by A) 

2. h(r) i A 

3. b + (r) \ A consists of exactly one element; we will denote it by a r . 

Our algorithm is based on the following result allowing us to restrict attention to the 
program P(A) (the statement of the lemma and its proof rely on the terminology 
introduced above). 

Lemma 4-3 

Let A be a set of atoms. A proper logic program P has an A-based stable model if 
and only if P(A) has an A-based stable model M = A U {a}, such that a £ {a r : r 6 
P'(A)}. 

Proof: (=>) Let M be an A-based stable model of P. Assume that M — A U {a}, 
for some a <£ A. Since P(A) M C P M , LM(P(A) M ) C LM(P M ) = M. Since M is 
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A-based, we have that M C LM{P(A) M ). It follows that M is an A-based stable 
model of P(A). 

Let us assume that there is a rule s G P'(A) such that a = a s . The rule s is not 
blocked by A. Since a G we have that a ^ 6 _ (s) (we recall that all rules in 

P are proper). Hence, s is not blocked by {a} either. Consequently, horn(s) G P M . 
Since s G P'(A), the body of horn(s) (that is, b + (sj) is contained in M. The set 
M is a least model of P M . In particular, M satisfies horn(s). Thus, it follows 
that h(s) G M. In the same time, h(s) ^ a (as s is proper). Thus, /i(s) G A, a 
contradiction (we recall that s G P'(A)). It follows that a ^ {a r :r G P'(A)}. 
(<=) We will now assume that M = A U {a} is an A-based stable model of P(A) 
such that a <£ {a r :r G P'(A)}. Similarly as before, we have M = LM(P(A) M ) C 
LM(P M ). Let us assume that LM(P M ) \ M ^ 0. Then there is a rule t in P M 
such that the body of t is contained in M and h(t) £ M. Let s be a rule in P 
that gives rise to t when constructing the reduct. Assume first that the body of t 
(that is, b + (s)) is contained in A. Then s G P (A), t G P(A) M and, consequently, 
h(t) G LM(P{A) M ) = M, a contradiction. 

Thus, the body of f is not contained in A. Since the body of t is contained in M, 
it consists of a and, possibly, some other elements, all of which are in A. It follows 
that s G P'{A). Consequently, a = a s and a G {a r :r G P'(A)}, a contradiction. 
Thus, LM(P M ) = M, that is, M is a stable model of P. Since M = LM(P(A) M ), 
it follows that M is an A-based model of P. □ 

Let A be a set of atoms. A logic program with negation, P, is an A-program if 
P = P(A), that is if for every rule rePwe have b+(P) C A and 6~(P) H A = 0. 
Clearly, the program P(A), described above, is an A-program. We will now focus 
on A-programs and their A-based stable models. 

Let A be a set of atoms. We denote by R(A) the set of all proper Horn rules 
over the set of atoms A. Clearly, the cardinality of R(A) depends on the cardinality 
of A only. Further, we define V(A) to be the set of all Horn programs Q C R{A) 
satisfying the condition LM(Q) — A. As in the case of R(A), the cardinality of 
V{A) also depends on the size of A only. 

We will now describe conditions that determine whether an A-program P has an 
A-based stable model. To this end, with every atom a G At(P) \ A, we associate 
the following values: 

• F(a) = 1 if there is a rule s in P with h(s) Au{a} and a £ b~(s); F(a) = 0, 
otherwise 

• G(a) = the number of rules s in P with h(s) — a and a £ b~(s). 

Further, with every proper Horn rule r G R(A) and every atom a G At{P) \ A, we 
associate the quantity: 

• H(r, a) = 1 if there is a rule s in P with horn(s) — r and a b~(s); H(r, a) = 
0, otherwise. 

The following lemma characterizes A-based stable models of an A-program. Both 
the statement of the lemma and its proof rely on the terminology introduced above. 

Lemma 4-4 
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Let A be a set of atoms, let P be an ^-program and let a be an atom such that 
a £ At(P)\A. Then 4u{a} is an A-based stable model of P if and only if F(a) = 0, 
G(a) > 0, and for some program Q £ V{A) and for every rule r £ Q, H(r, a) > 0. 

Proof: (=>) We denote M — A U {a} and assume that M is an A-based stable 
model for P. It follows that M = LM(P M ). Let P A be the subpro gram of P 
consisting of those rules of P whose head belongs to A. Since M is an A-based 
stable model of P, we have A — LM{P A vI ). Let Q be the program obtained from 
by removing multiple occurrences of rules. Clearly, Q £ V(A). It follows directly 
from the definition of the reduct that for every rule r £ Q, H(r, a) = 1. 

Next, we observe that a £ LM(P M ). Thus, G(a) > 0. Let us assume that F(a) = 
1. Let r be a rule in P such that h(r) A U {a} and a b~(r). Since P is an A- 
program, A n b~(r) = 0. Thus, it follows that horn(r) £ P M . We also have that 
b + (r) CiC M. Since M is a model of P M , /i(r) £ M . However, in the same time 
we have that h{r) £ A U {a}(= M), a contradiction. It follows that F(a) = 0. 
(^) We now assume that for some a £ At(P) \ A, F(a) — 0, G(a) > and for 
some program Q £ V(A) and for every rule r £ Q, H(r, a) = 1. As before, we set 
M = A U {a}. We will show that M = LM{P M ). 

First, since P is an A-program and H(r, a) = 1 for every rule r £ Q, it follows 
that Q C P(A) M . Thus, A C LM(P(A) M ). Second, we have that G(a) > 0. Thus, 
there is a rule r £ P such that /i(r) = a and a ^ b~{r). It follows that horn(r) £ Q 
and /iorn(r) G P M . Since Q C P{A) M , A = LM(Q) and 6+(r) C A, we obtain 
that a G LM(P(A) M ). Thus, M C L¥(P(i) M ). Finally, since F(a) = 0, we have 
that for every rule s £ P such that a g b~(s), h(s) £ M. Thus, LM(P M ) does not 
contain any atom not in M. Consequently, M = LM(P M ) and M is a stable model 
of P. Since M C LM(P(A) M ), M is an A-based stable model of P. □ 

We will discuss now effective ways to compute values F(a), G{a) and H(r,a). 
Clearly, computing the values G(a) can be accomplished in time linear in the size 
of the program, that is, in time 0(m). Indeed, we start by initializing all values 
G(a) to 0. Then, for each rule s £ P, we set G(h(s)) := G(h(s)) + 1 if h(s) <£ b~(s), 
and leave G(h(s)) unchanged, otherwise. To decide which is the case requires that 
we scan all negated lierals in the body of s. That takes time 0(\b~(s)\). Thus, the 
overall time is 0(m). 

Computing values F(a) and H(r,a) is more complicated. First, we prove the 
following lemma. 

Lemma 4-5 

Let P be an A-program, let a £ At(P) \ A and let r £ R(A). Then 

1. F(a) = 1 if and only if a <£ f){{h(s)} U b~(s) : s £ P,h(s) £ A}. 

2. H(r, a) = 1 if and only if a £ f){b~( s ) ■ s £ P, horn(s) = r}. 

Proof: (1) Let us assume first that F(a) = 1. Then there is a rule s £ P such 
that h(s) £ All {a} and a £ b~(s). Thus, a £ {h(s)} U b~(s). Consequently, the 
identity a £ n{{M s )} u b~(s) : s £ P,h(s) £ A} follows. All the implications in 
this argument can be reversed. Hence, we obtain the assertion (1). 
(2) Let us assume that H(r, a) = 1. Then, there is a rule s £ P such that horn(s) = r 
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and a £ b~(s). Consequently, a ^ n{^( s ) : s £ P,horn(s) = r}. As in (1), all the 
implications are in fact equivalences and the assertion (2) follows. □ 
Lemma 4.5 shows that to compute all the values F(a) one has to compute the 
set 

f]{{h(s)}Ub-(s):s€P,h(s)^A}. 

To this end, for each atom a we will compute the number of sets in {{h(s)}Ub~ (s) : 
s G P, h(s) ^ A} that a is a member of. We will denote this number by C(a). We 
first initialize all values C(a) to 0. Then, we consider all sets in {{h(s)} U b~(s) : 
s G P, h(s) f. A} in turn. For each such set and for each atom a in this set we 
set C{a) := C{a) + 1. The set f]{{h(s)} U b~(s) : s G P, h(s) £ A} is given by all 
those atoms a for which C(a) is equal to the number of sets in {{h(s)} U b~(s) : 
s G P, h(s) f. A}. It is clear that the time needed for this computation is linear in 
the size of the program (assuming appropriate linked- list representation of rules). 
Thus, all the values F (a) can be computed in time linear in the size of the program, 
that is, in O(m) steps. 

To compute values H(r,a) we proceed similarly. First, we compute all the sets 
{s : s G P,horn(s) = r}, where r G R(A). To this end, we scan all rules in P in 
order and for each of them we find the rule r G R(A) such that horn(s) = r. Then 
we include s in the set {s : s G P, horn(s) = r}. Given s, it takes 0(.g|A|) steps to 
identify rule r (where g is some function). Indeed, the size of b + {s) is bound by \A\ 
as P is an yl-program. Moreover, the number of rules in R(A) depends on \A\ only. 
Thus, the task of computing all sets {s : s G P, horn(s) = r}, for r G R(A), can be 
accomplished in 0(g(|A|)|P|) steps. Next, for each these sets of rules, we proceed as 
in the case of values F(a), to compute their intersections. Each such computation 
takes time O(m), where m = size(P)). Thus, computing all the values H (r, a) can 
be accomplished in time 0(g(|A|)|P| + \R(A)\m) = 0(f(\A\)m), for some function 
/• 

We can now put all the pieces together. As a result of our considerations, we 
obtain the following algorithm for deciding the problem SSM(k). 



Algorithm to decide the problem SSM(k), k > 1 
Input: A logic program P (k is not a part of input) 

(0) if is a stable model of P then return YES and exit; 

(1) P := the set of proper rules in P; 

(2) for every A C At(P) with \A\ < k - 1 do 

(3) compute the set of rules R(A) and the set of programs V(A); 

(4) compute the program P(A); 

(5) compute the program P'(A) and the set B = {a r : r G P'(A)}; 

(6) given P(A) and R(A), compute tables F, G and H (as described above); 

(7) for every a G At(P(A)) \ A \ B do 

(8) if 

(9) F(a) = 0, G{a) > and 

(10) there is a program Q G V(A) s. t. for every rule r G Q, H(r, a) > 
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(11) then report YES and exit; 

(12) report NO and exit. 



The correctness of this algorithm follows from Lemmas 4.2 - 4.4. We will now 
analyze the running time of this algorithm. Clearly, line (0) can be executed in 
O(m) steps. As we already observed, rules that are not proper can be eliminated 
from P in time 0(m). Next, there are 0(n k ~ 1 ) iterations of loop (2). In each of 
them, line (3) takes time 0(/i(fc)), for some function /i (let us recall that |-R(A)| 
and ('P(A)I depend on \A\ only). Further, lines (4) and (5) can be executed in 
time 0{m). Line (6), as we discussed earlier, can be implemented so that to run in 
0(f(k)m) steps. Loop (7) is executed 0(n) times and each iteration takes 0(/a(A;)) 
steps, for some function /2 (let us again recall that ^(A) | depends on k only). 
Thus, the running time of the whole algorithm is 0(F(k)mn k ~ 1 ), for some integer 
function F. Consequently, we get the following result. 

Theorem 4-6 

There is an integer function F and an algorithm A such that A decides the problem 
SSM(k) and runs in time 0(F(k)mn ) (the constant hidden in the "big Oh" 
notation does not depend on k). 



5 Complexity of the problem SSM 

The algorithm outlined in the previous section is not quite satisfactory. Its run- 
ning time is still high. A natural question to ask is: are there significantly better 
algorithms for the problems SSM(k)? In this section we address this question by 
studying the complexity of the problem SSM. Our goal is to show that the problem 
is difficult in the sense of the W hierarchy. We will show that the problem SSM is 
M / [2]-hard and that it is in the class W[3]. To this end, we define the (< k) -weighted 
t-normalized satisfiability problem as: 

WS-{t) Given a i-normalizcd formula ip, decide whether there is a model of tp with 
at most k atoms (k is a parameter). 

The problem WS-(t) is a slight variation of the problem WS(t). It is known to be 



complete for the class W[t], for t > 2 (see ( Downey fc Fellows, 1997 ), page 468). To 
show W[2]-hardness of SSM, we will reduce the problem WS-(2) to the problem 
SSM. Given the overwhelming evidence of fixed-parameter intractability of prob- 



lems that are VF[2]-hard (Downey & Fellows, 1997), it is unlikely that algorithms for 



problems SSM(k) exist whose asymptotic behavior would be given by a polynomial 
of order independent of k. To better delineate the location of the problem SSM in 
the W hierarchy we also provide an upper bound on its hardness by showing that 
it can be reduced to the problem WS— (3), thus proving that the problem SSM 
belongs to the class W[3]. 

We will start by showing that the problem SSM(k) is reducible (in the sense of 
the definition from Section ||) to the problem WS-(3). To this end, we describe 
an encoding of a logic program P by means of a collection of clauses T(P) so that 
P has a stable model of size at most k if and only if T(P) has a model with no 
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more than (k + l)(fc 2 + 2k) atoms. In the general setting of the class NP, an explicit 
encoding of the problem of existence of stable models in terms of propositional 
satisfiability was described in (Bcn-Eliyahu & Dcchtcr, 1994). Our encoding, while 
different in key details, uses some ideas from that paper. 

Let us consider an integer k and a logic program P. For each atom q in P 
let us introduce new atoms c(q), c(q,i), 1 < i < k + 1, and c~(q,i), 2 < i < 
k + 1. Intuitively, atom c(q) represents the fact that in the process of computing 
the least model of the reduct of P with respect to some set of atoms, atom q is 
computed no later than during the iteration k + 1 of the van Emden-Kowalski 
operator. Similarly, atom c(q, i) represents the fact that in the same process atom q 
is computed exactly in the iteration i of the van Emden-Kowalski operator. Finally, 
atom c~ {q, i), expresses the fact that q is computed before the iteration i of the van 
Emden-Kowalski operator. The formulas Fi(q,i), 2 < i < fc + 1, and ^2(9) describe 
some basic relationships between atoms c(q), c(q, i) and c~ (q, i) that we will require 
to hold: 

Fi (q, i) = c~ (q, i) ^ c(q, 1) V ... V c(q, 

F 2 {q) = c{q) & c(q, 1) V ... V c(q, k + 1). 
Let r be a rule in P with h(r) — q, say 

r = q <— 01, . . . , a s , not(6i), . . . , not(6 t ). 

We define a formula Fs(r, i), 2 < i < k + 1, by 

F 3 (r,i) = c^(a 1 ,i) A ... A c~(a s ,i) A ->c(fei) A ... A ~^c(b t ) A -*c~(q,i). 

We define Fs(r, 1) = false (false is a distinguished contradictory formula in our 
propositional language) if s > 1. Otherwise, wc define 

F 3 (r,l) = -.c(6i)A...A-.c(6t). 

Speaking informally, formula F^(r, i) asserts that q is computed by means of rule r 
in the iteration i of the least model computation process and that it has not been 
computed earlier. 

Let ri, . . . ,r v be all rules in P with atom q in the head. We define a formula 
F±{q,i), 1 < i < k + 1, by 

%j)=c( g ,i)«f 3 (r 1 ,i)V...V% 1 i). 

Intuitively, the formula F^q, i) asserts that when computing the least model of the 
reduct of P, atom q is first computed in the iteration i. 

We now define the theory Tq(P) that encodes the problem of existence of small 
stable models: 

T (P) = {F 1 (q,i):qGAt(P),2<i<k + l}U{F 2 (q):qeAt(P)}U 
{F 4 (<7, i): q G At(P), 1 < i < k + 1}. 

Next, we establish some useful properties of the theory Tq(P). First, we consider 
a set U of atoms that is a model of T (P) and define 

M(U) = {qe At(P):c(q) £ U}. 
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Lemma 5.1 

Let U be a model of Tq(P) and let q G M(U). Then there is a unique integer i, 
1 < i < k + 1, such that c(g, i) G J7. 

Proof: Since {/ is a model of a formula -^2(9), there is an integer i, l<i<fc+l, 
such that c(<7, i) G {/. To prove uniqueness of such i, assume that there are two 
integers j\ and j 2 , 1 < ji < 3i < & + 1, such that c(q,ji) G £/ and c(<7, j 2 ) G ?7. 
Since £/ |= i*4(g, J2), it follows that there is a rule r G P with /i(r) = <? and such that 
U \= Fs(r,j2)- In particular, U (= ->c _ (g,j2). In the same time, since c(q,ji) G J7 
and [/ |= Fi(q,j2), we have c~(q,j2) G J7, a contradiction. □ 
For every atom g G M(U) define z g to be the integer whose existence and unique- 
ness is guaranteed by Lemma 5T. Define ijj = max{i g : q G Af ([/)}. Next, for each 



i, 1 < i < i(7, define 

[M(Cr)] i = {g€Af(lO:i g = i}. 

Lemma 5.2 

Let £7 be a model of Tq(P). Under the terminology introduced above, for every i, 
l<i<iu, [M(U)]i + 0. 

Proof: We will proceed by downward induction. By the definition of iu, [M {U)\i u ^ 
0. Consider i, 2 < i < iu, and assume that [M (£/)], ^ 0. We will show that 
[M(U)] l -i ^ 0. Let q G [M(U)]i. Clearly, c(q,i) G U and, since U \= F A (q,i), 
there is a rule r = q <— ai, . . . , a s , not(&i), . . . , not(6t) such that U \= F 3 (r 7 i). 
Consequently, for every j, 1 < j < s, c — (aj,i) G £7. Assume that for every j, 1 < 
j < s, c~(a,j, i — 1) G £7. Since [/ |= c~ (q, i — 1) c~((j, i) and since £7 |= -<c~(q, i), 
it follows that £7 \= ~^c~ (q, i — 1). Consequently, £7 satisfies the formula Fs(r, i — 1) 
and, so, £7 |= F4{q,i — 1). It follows that c(q,i — 1) G £7, a contradiction (we recall 
that i 9 = i). Hence, there is j, 1 < j < s, such that c(a,j,i — 1) G £7. It follows that 
Oj G [M(Z7)]i_i and [M(Z7)] i _ 1 7^ 0. □ 

Lemma 5.3 

Let [7 be a model of T (P) and let \M(U)\ < k. Then 

1. ijj < k, and 

2. M(E7) is a stable model of P. 

Proof: (1) The assertion follows directly from the fact that \M(U)\ < k and from 
Lemma |5.2j . 

(2) We need to show that M(£7) = LM{P M W). We will first show that M(E7) C 
LM{P M W). Since M(£7) = \JZii M (U)h, we will show that for every i, 1 < i < 
iu, [M(U)]i C LM{P M W). We will proceed by induction. Let q G [M(C7)]i. It 
follows that there is a rule r such that £7 |= ^(r, 1). Consequently, r is of the 
form r = q <— not (61), . . . , not(& t ) and £7 |= _, c(6i) A ... A ^c(b t ). Hence, for every 
j, 1 < j < t, bj £ M(U). Consequently, the rule (q <— .) is in p M ( u ) and, so, 
q G LM(P M W). The inductive step is based on a similar argument. It relies on 
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the inequality ijj < k we proved in (1). We leave the details of the inductive step 
to the reader. 

We will next show that LM(P M W) C M{U). We will use the characterization of 
LM(P M W) as the limit of the sequence of iterations of the van Emden-Kowalski 
operator T p m(u) : 

oo 

LM(P M ^) = {jT pM ^($). 

i=0 

We will first show that for every integer i, < i < k+1, we have: T pM{u) (0) C M (U) 
and for every q G T pU(u) (0), i q < i. 

Clearly T pM(u) (0) = C M(U). Hence, the basis for the induction is established. 
Assume that for some i, < i < k, Tp M([J) (0) C M(U) and that for every q G 
T pM{ u } (9), i q < i- Consider q G T l p \} {u) {%). If U \= c~(q,i+ 1), then c{q,v) G U for 
some v, 1 < v < i. Since U \= ^2(9), c(q) G U and q G M(U). By Lemma [Tl], it 
follows that i q = v. Hence, i q < i+ 1. 

Thus, assume that [7 |= ^c~(q,i + 1). Since q G T p + ^ {u) (%), there is a rule 

r = 9 <— ai, . . . , a s , not(&i), . . . , not(6 4 ) 

in P such that bj M (U), for every j, 1 < j < t, and a,j G T pM{u) (0), 1 < i < s. By 
the induction hypothesis, for every j, 1 < j < s, we have a, G M(U) and i aj < i. 
It follows that [/ (= Fs(r, i + 1) and, consequently, that c(q,i + 1) G U. Since 
[7 |= F 2 (q), c(q) G [7 and q G M(t7). It also follows (Lemma that = i + 1. 

Thus, we proved that U;=o T pM(t/) (0) C Af([7). Since |M(E7)| < fc, there is j, 
< j < k such that T pM{u) {%) = T p \} {u) (0). It follows that for every j', j < f, 

T PM{U) (9) = T J pM{u) {%). Consequently, T pM(u) (9) C A7([7) for every non-negative 
integer i. □ 
Consider now a stable model M of the program P and assume that \M\ < k. 
Clearly, M = \J°° =1 T pM {%). For each atom q G Af define s q to be the least integer 
s such that q G T pM {%). Clearly, s q > 1. Moreover, since \M\ < k, it follows that 
for each q G M, s q < k. Now, define 

U M = {c(q),c(q,s q ):q G A/} U {<T(>,i):q G M, s q < i < k + 1} 
Lemma 5.4 

Let M be a stable model of a logic program P such that \M\ < k. Under the 
terminology introduced above, the set of atoms Um is a model of Tq(P). 

Proof: Clearly, U M |= Pi (q,i) for q G ,4i(P) and 2 < i < k + 1, and [7 M |= P 2 (q) 
for q G ,4i(P). 

We will now show that Um h ^ife *)> f° r 1 e -^(-P) an d i = 1, 2, . . . , First, 
we will consider the case q G M. There are three subcases here depending on the 
value of i. 

We start with i such that s q < i < k + 1. Then Um ~' c_ ( ( Zi*)- ^ follows 
that Um \£ -^3(7, *) for every rule r G P such that fc(r) = q. Since Um c (9;*); 
(7m 1=^(9,0. 
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Next, we assume that i — s q . Then, there is a rule r = q <— a±, . . . , a s , not (61), . . . , 
not(bf) in P such that bj £ M, for every j, 1 < j < t, and aj G T^}^), 1 < j < s. 
Clearly, Um \= F 3 (r, i). Since Um \= c(q, i), it follows that Um H Fi(q, f° r * = -V 

Finally, let us consider the case 1 < i < s q . Assume that there is rule r G P such 
that h(r) = q and Um \= F^(r, «)■ Let us assume that r = q *— ax, . . . , a s , not (61), . . . , 
not(6 t ). It follows that for every j, 1 < j < t, Um \= _, c(6 J ). Consequently, for every 
j, 1 < j < t, bj ^ M and the rule r' = q <— 01, . . . , a s belongs to the reduct P M . In 
addition, for every j , 1 < j < s, c~(aj,i) 6 Um- Thus, aj G M and s ai < i — 1 . This 
latter property is equivalent to aj S T l ~ h }{%). Thus, it follows that (7 G Tp M (0) and 
s q < i — a contradiction with the assumption that i < s q . Hence, for every rule r 
with the head q, Um F3(r,i)- Since for i < s q , c(q,i) Um, Um \= F±(q,i). 

To complete the proof, we still need to consider the case q £ M . Clearly, for 
every i, 1 < i < k + 1, Um V= c(q,i). Assume that there is i, 1 < i < k + 1, 
and a rule r such that h(r) — q and Um \= F3{r,i). Let us assume that r is of 
the form q *— a\, . . . , a s , not (61), . . . , not (64). It follows that c~(aj,i) G C/jif and, 
consequently, a 3 G Af for every j, 1 < J < s. In addition, it follows that for every j, 
1 < j < t, Um \= ->c(bj) and, consequently, bj M. Thus, q <— a\, . . . , a s belongs 
to the reduct P M and, since M is a model of the reduct, q G M, a contradiction. 
It follows that for every i, 1 < i < k + 1, Um \= F±{q, i). □ 

For each atom q G At(P), let us introduce k 2 + 2k new atoms d(q,i), 1 < i < 
k 2 + 2k, and define 

T{P) = T (P) U {c(g) ^ d(g, i): 1 < i < k 2 + 2k}. 



Lemmas 5.1 - 5.4 add up to a proof of the following result. 



Theorem 5.5 

Let k be a non-negative integer and let P be a logic program. The program P has 
a stable model of size at most k if and only if the theory T(P) has a model of size 
at most (k + \){k 2 + 2k). 



Proof: (=£-) Let M be a stable model of P such that |M| < k. By Lemma 5.4, the 
set Um is a model of Tq(P) Consequently, the set 

U = U M U {d(g, i): q G M, 1 < i < fc 2 + 2fc} 

is a model of T(P). Moreover, it is easy to see that \Um\ < 2k + k 2 . Hence, \U\ < 
2k + k 2 + k{k 2 + 2k) = (k + l){k 2 + 2k). 

Conversely, let us assume that some set V, consisting of atoms appearing in T(P) 
and such that \V\ < (k + l)(k 2 + 2k), is a model of T(P). Let us define U to consist 
of all atoms of the form c(q), c(q,i) and c~(q,i) that appear in V. Clearly, U is 
a model of Tq(P). Let us assume that M(U) > k + 1 (we recall that the notation 



M(U) was introduced just before Lemma 5.1 was stated). Then, there are at least 



(k + l)(k 2 + 2k) atoms of type d(q, i) in V. Consequently, V > (k + l)(k 2 + 2k) as 
it contains also at least k + 1 atoms c(q), where q G M{U). This is a contradiction. 
Thus, it follows that \M(U)\ < k. Moreover, by Lemma [53| , M(U) is a stable model 
of P. □ 
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Let us now define the following sets of formulas. First, for each atom q G At(P) 
we define 

C (q) = {-<c(q) V i): 1 < i < k 2 + 2k} U {c(g) V -.<%, i): 1 < i < k 2 + 2k}. 
Next, we define 

d(q,i) = {^c~(q, i) V c(q, 1) V ... V c(q, i — 1)} U {->c(q,j) V c~(q, i): I < j <i — 1}, 

C 2 (q) = {-c(g) V c(q, 1) V ... V c(g, fc + 1)} U {^c(q, j) V c(q): l<j<k+l}, 
and 

C 4 (g,i) = {-■<:(?, i) VF 3 (n,i) V ... VF 3 (r„i)}U {->F 3 (rj, i) V c(g, i): 1 < j < «}, 

where {ri, . . . , r^} is the set of all rules in P with q in the head. 
Clearly, the theory 

T C (P) = {C (g):«e At(P)}U{(7i(g,i):geAt(P),2<i<fc+l}U 
{C 2 (g): q G At(P)} U {C 4 (q, i): q G At(P), 1 < i < fc + 1} 

is equivalent to the theory T(P). Moreover, it is a collection of sums of products of 
literals. Therefore, it is a 3-normalized formula. By Theorem it follows that the 
problem 55Mcan be reduced to the problem WS-(3). Thus, we get the following 
result. 

Theorem 5.6 

The problem SSM(k) G W[3\. 

Next, we will show that the problem WS-(2) can be reduced to the problem 
SSM. Let C = {ci, . . . , Cp} be a collection of clauses. Let A — {xi, ... ,x r } be the 
set of atoms appearing in clauses in C . For each atom x £ A, introduce k new 
atoms x(i), 1 < i < k. By Si, 1 < i < k, we denote the logic program consisting of 
the following n clauses: 

x\(i) <— not(a;2(i)), . . . , not(x r {i)) 

x r (i) <— not(.xi (i)), . . . , not(av_i(i)) 

Define S = UiLi Clearly, each stable model of S is of the form {xj t (1), . . . , xj k (k)}, 
where 1 < j p < r for p = 1, . . . , k. Sets of this form can be viewed as representations 
of nonempty subsets of the set A that have no more than k elements. This repre- 
sentation is not one-to-one, that is, some subsets have multiple representations. 
Next, define Pi to be the program consisting of the clauses 

xj < Xj(i), j = 1, . . . ,r, i = 1, 2, . . . , k. 

Stable models of the program S U Pi are of the form {^(l), . . . ,Xj k (k)} U M, 
where M is a nonempty subset of A such that \M\ < k and Xj t , . . . , Xj k enumerate 
(possibly with repetitions) all elements of M. 
Finally, for each clause 

c = a x V . . . V a s V ->&i V . . . V -^b t 
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from C define a logic program clause p{c) : 

p{c) = f <- h, ■ ■ ■ , h, not(ai), . . . , not(a s ), not(/) 
where / is yet another new atom. Define P 2 = {p(c):c & C} and P c = S 1 U Pi U P 2 . 
Theorem 5. 7 

A set of clauses C has a nonempty model with no more than k elements if and only 
if the program P c has a stable model with no more than 2k elements. 

Proof: Let M be a nonempty model of C such that \M\ < k. Let Xj 1 , . . . , Xj k 
be an enumeration of all elements of M (possibly with repetitions). Then the set 
M' = {xj 1 (1), . . . , Xj k (k)} U M is a stable model of the program S U Pi. Since M 
is a model of C, it follows that (P C ) M ' = (S U Pi) M ' U F, where F consists of the 
clauses of the form 

/ <-&i,.. .,b t , 

such that t > 1 and for some j, 1 < j < t, bj £ M' . Since M' = LM((S U Pi) M '), 
it follows that 

W = LM({S U Pi) M ' U P) = LM({P C ) M '). 

Thus, M' is a stable model of P c . Since |M'| < 2k, the "only if" part of the 
assertion follows. 

Conversely, assume that M' is a stable model of P c . Clearly, / ^ M' . Conse- 
quently, 

LM((5UPi) M ') = LM((5U?i UP 2 ) M ') = LM((P C ) M ') = M'. 

That is, M' is a stable model of S U Pi. As mentioned earlier, it follows that 
M' = {xj 1 (1), . . . , Xj k (k)} U M, where M is a nonempty subset of At(P) such that 
\M\ < fc and Xjj, . . . , a:^ is an enumeration of all elements of M. 

Consider a clause c = a\ V . . . V a s V — 161 V ... V from C. Since M' is a stable 
model of P c , it is a model of P c . In particular, M' is a model of p(c). Since / ^ M', 
it follows that M' \= c and, consequently, M [= c. Hence, M is a model of C. □ 

Now the reducibility of the problem WS-{2) to the problem SSM is evident. 
Given a collection of clauses C, to check whether it has a model of size at most 
k, we first check whether the empty set of atoms is a model of C. If so, we return 
the answer YES and terminate the algorithm. Otherwise, we construct the program 
P c and check whether it has a stable model of size at most 2k. Consequently, we 
obtain the following result. 

Theorem 5.8 

The problem SSM is W[2]-hard. 

6 Open problems and conclusions 

The paper established several results pertaining to the problem of computing small 
and large stable models. It also brings up interesting research questions. 

First, we proved that the problem LSM is in the class FPT. For problems that 
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are fixed-parameter tractable, it is often possible to design an algorithm running 
in time 0{p(N) + f(k)), where N is the size of the problem, k is a parameter, p 
is a polynomial and / is a function (Downey & Fellows, 1997). Such algorithms 
are often practical for quite large ranges of N and k. The algorithm for the LSM 
problem presented in this paper runs in time 0(m2 k+k ). It seems plausible it can 
be improved to run in time 0(m + f(k)), for some function /. Such an algorithm 
would most certainly be practical for wide range of values of m and k. We propose 
as an open problem the challenge of designing an algorithm for computing large 
stable models with this time complexity. 

There is a natural variation on the problem of computing large stable models: 
given a logic program P and an integer k (parameter), decide whether P has a 
stable model of size at least |j4t(P)| — k. This version of the problem LSM was 
recently proved by Zbigniew Lone and the author to be W[3]-hard (and, hence, 
fixed-parameter intractable) ( |Lonc fc Truszczynski, 2000 ). The upper bound for 
the complexity of this problem remains unknown. 

In the paper, we described an algorithm that for every fixed k, decides the exis- 
tence of stable models of size at most k in time 0(n k ~ 1 m), where n is the number of 
atoms in the program and m is its size. This algorithm offers only a slight improve- 
ment over the straightforward "guess-and-check" algorithm. An interesting and, it 
seems, difficult problem is to significantly improve on this algorithm by lowering 
the exponent in the complexity estimate to ah, for some constant a < 1. 

We also studied the complexity of the problem SSM and showed that it is fixed- 
parameter intractable. Our results show that SSM is Vy[2]-hard. This result implies 
that the problem SSM is at least as hard as the problem to determine whether 
a CNF theory has a model of cardinality at most k, and strongly suggests that 
algorithms do not exist that would decide problems SSM(k) and run in time 0(n c ), 
where c is a constant independent on k. For the upper bound, we proved in this 
paper that the problem SSM belongs to class W[3]. Recently, Zbigniew Lone and 
the author (Lone & Truszczynski, 200C) showed that the problem SSM is, in fact, 
in the class W[2]. 
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